gusucode.com > 现代通信系统——使用MATLAB（+全部程序） > 现代通信系统——使用MATLAB（+全部程序）/现代通信系统——使用MATLAB（+全部程序）/Matlab 程序/Chapter5/ip_05_06.m

    % MATLAB script for Illustrative Problem 5.6. 
echo on
SNRindB1=0:1:15;
SNRindB2=0:0.1:15;
for i=1:length(SNRindB1),
   smld_err_prb(i)=smldPe56(SNRindB1(i));   	% simulated error rate
   echo off;
end;
echo on;
for i=1:length(SNRindB2),
   SNR=exp(SNRindB2(i)*log(10)/10);         	% signal-to-noise ratio
   theo_err_prb(i)=Qfunct(sqrt(SNR/2));  	% theoretical error rate
   echo off;
end;
echo on;
% Plotting commands follow.
semilogy(SNRindB1,smld_err_prb,'*');
hold
semilogy(SNRindB2,theo_err_prb);